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We propose and study a simple model of dynamical redistribution of capital in a diversified portfolio. 
We consider a hypothetical situation of a portfolio composed of N uncorrelated stocks. Each stock 
price follows a multiplicative random walk with identical drift and dispersion. The rules of our model 
naturally give rise to power law tails in the distribution of capital fractions invested in different 
stocks. The exponent of this scale free distribution is calculated in both discrete and continuous 
time formalism. It is demonstrated that the dynamical redistribution strategy results in a larger 
typical growth rate of the capital than a static "buy-and-hold" strategy. In the large N limit the 
typical growth rate is shown to asymptotically approach that of the expectation value of the stock 
price . The finite dimensional variant of the model is shown to describe the partition function of 
directed polymers in random media. 



I. INTRODUCTION 

The problem of finding an investment strategy with the 
best long-term growth rate of the capital is of tremendous 
practical importance. The traditional theory of portfo- 
lio optimization is stationary in origin It answers 
the question of optimal distribution of the capital be- 
tween different assets (optimal asset allocation), but in 
general gives no prescription on how to maintain this op- 
timal allocation at all times. In this work we propose a 
simple model of dynamical allocation of capital. Some- 
what countcrintuitivcly, in order to optimize the growth 
rate of the capital an investor has to sell assets which 
have increased in price since the last update, and buy 
those which have decreased. In doing so he sells stocks 
when they are "overpriced" and buys them when they 
are "underpriced" , which is clearly advantageous. As we 
demonstrate below, in our model an investor who actively 
manages his portfolio in such a fashion almost certainly 
does better than one who follows a static "buy-and-hold" 
strategy. 

The nontrivial properties of the problem come 
from the multiplicative nature of stock price fluctua- 
tions. Throughout this manuscript we assume that on 
timescales of interest to us the prices of individual assets 
follow a multiplicative random walk. In other words, the 
ratio of stock prices at two consecutive times, at which 
the investor buys or sells stock, is a random number, un- 
correlated with the current price and with the history of 
price changes in the past. There are many peculiarities of 
such noisy multiplicative dynamics, especially regarding 
expectation values of random variables. Traditional ex- 
pectation (average) value is of little relevance here. The 
reason for this is that the dominant contribution to the 
expectation value of a random variable subject to mul- 
tiplicative noise comes from exponentially unlikely out- 
comes when the variable is exponentially large. For any 
finite number of realizations (and in real world one always 
deals with just one realization) this expectation (average) 



value is very unlikely to appear. On the other hand, the 
typical value of such random variable, defined as the me- 
dian of its probability distribution, constitutes a more 
realistic property. 

Just like in the static portfolio theory, our strategy fa- 
vors the diversification, i.e. increasing the number of as- 
sets in the portfolio. We demonstrate that in our model 
the diversification reduces fluctuations, and makes the 
growth rate of the typical value of the capital to be closer 
to that of its expectation value. However, for any finite 
number of assets, these two growth rates are still differ- 
ent. 

Under the rules of dynamical redistribution of funds, 
which we employ in this manuscript, the distribution of 
shares of the total capital invested in individual assets 
naturally acquires a power law tail. This adds yet an- 
other example of how a scale free distribution can arise 
out of multiplicative dynamics without fine-tuning of any 
sort. We derive the analytical expression for the expo- 
nent r of this power law. Somewhat surprisingly, in the 
weak coupling limit, corresponding to slow redistribution 
of funds between the assets, this exponent has a "supe- 
runiversal" value r = 2. It gradually increases with the 
coupling constant and becomes infinite in the limit where 
the capital is equally redistributed between assets after 
each time step. 

The rules of redistribution of capital can be interpreted 
as fully-connected (infinite-dimensional) limit of the well 
known statistical model of directed polymers in the pres- 
ence of quenched disorder. This provides a new and excit- 
ing link between the physics of finance, and the problems 
lying on the forefront of modern theoretical condensed 
matter physics. 

The plan of the manuscript is as follows: to streamline 
the following introduction of our basic model, in Section 
II we review the well known (and not so well known) 
properties of a stochastic multiplicative dynamics. We 
remind the reader the formulas for average and typical 
value of a single multiplicative random walk, formulate 
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the "continuous" time approach to this problem, and re- 
fresh in reader's memory the formalism of Ito stochas- 
tic calculus, necessary for our purposes. Then we re- 
view recent results on natural appearance of power law 
distributions in a situation when a single multiplicative 
random walk is pushed against lower wall [^|j3| , prevent- 
ing the random variable from falling below certain value. 
Finally, we describe the multiplicative stock price and 
capital dynamics used throughout this manuscript. 

In Section III we analyze the behavior of the typ- 
ical and average values of the capital in a "buy-and- 
hold" strategy, where the capital was initially equally 
distributed between N independent assets with the same 
typical growth rate and dispersion, and no further re- 
distribution ever took place. We demonstrate that after 
a logarithmically short initial period of time, the typical 
growth rate of the capital is limited to the typical growth 
rate of the price of the assets, and is significantly smaller 
than their average (expectation) growth rate (or average 
return per capital of this asset) 

In Section IV we show that the growth rate of investor's 
capital can be significantly increased by following an ac- 
tive, dynamic redistribution strategy. In this strategy 
at each time step the investor sells some shares of every 
stock with current value of invested capital above the all- 
stock average, and buys some shares of every stock below 
this average. We analyze the consequences of this strat- 
egy in both discrete and continuous time formalisms and 
demonstrate that in both cases these rules naturally give 
rise to a scale free distribution of fractions of individ- 
ual stock capitals in the total capital. We proceed with 
deriving analytical expressions for the exponent of this 
distribution, and the typical growth rate of the capital 
in this situation. This rate for our strategy proves to be 
larger than that in the static "buy-and-hold" strategy, 
which a posteriori justifies our approach. However, as 
should be expected, the total capital is still subject to 
the multiplicative noise, and therefore its typical growth 
rate is still smaller than the average growth rate. We 
demonstrate that in the limit N — > oo these two rates 
asymptotically converge as some power of 1/JV. 



II. REVIEW OF RESULTS FOR A SINGLE 
MULTIPLICATIVE RANDOM WALK 

A. Typical and average values of a multiplicative 
random walk 

Consider a stochastic process in which at each time 
step a variable W(t) is multiplied by a positive random 
number e* 7 ^ , where 77(f) is drawn from some probability 
distribution ^(77): 



W(f + 1) = e v ® W(t). 



(1) 



We adopt the initial condition W(Q) = 1. For the new 
variable h(t) = lnVF(f) this process is just a random 



walk with an average drift v = (77) and a dispersion D = 
(rj 2 ) — (77) 2 . The corresponding equation of motion is 
simply 



h(t + 1) = h(t)+ 77(f). 



(2) 



In recent literature it has been observed that average 
and typical values of W(t) in such a process can be very 
different. One of the precise definitions of the typical 
value of a random variable is the median [f| of its prob- 
ability distribution, i.e. for Wt yp one has the property 
that Prob(PF > W typ ) = Prob(IT < W typ ) = 1/2. By 
definition W typ (t) = e 1 ""!'). 

The central limit theorem implies that asymptotically 
the distribution P{h, t) can be approximated with a 
Gaussian 
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Therefore, the median (as well as average and most 
probable values) of h(t) changes linearly with time, and 
the rate of this change is given by the drift velocity 
v = (77) of the corresponding random walk: In Wt yp (f) = 
(lnW(t)) = {r,)t. 

On the other hand the expectation (average) value of 
W{t) changes as (W(t + 1)) = (e"i){W(t)) (since 77(f) and 
W(t) are uncorrelated) . Hence, ln(lT(f)) = ln(e r? ) f also 
depends linearly on time but with a different slope. It 
is easy to show that for any distribution ln(e'') > (77), 
so that the average value of W always grows faster than 
its typical value and after some time one has (W(t)) 
Wtyp(f). This exponentially large discrepancy between 
typical and average values of W is due to the long tails 
of P(W,t), but the events constituting these tails are 
extremely rare. 

For future use we derive analytic expressions for the 
growth rate of (W m (t)) in a simple case, when 77 is 
drawn from a Gaussian distribution with average value 
v = (77) and dispersion D = (rj 2 ) — (77) 2 . Since the dy- 
namics of W m is given by W m (t + 1) = e^ m W m (t), 
for (W m (t)) one has (W m (t)) = (e" ro )*. The integral 
JZo dr J eVm e-^- v ) 2 / 2D /y/^dD can be easily taken and 
is equal to e m ( v + Dm / 2 ) . Therefore, for a Gaussian distri- 
bution one has 



(e' ?m ) 1 ^ m = e ' u +- Dm /2 
(W m {t)) = e m ( v + Dm / 2 )t, 



(4) 
(5) 

It is important to mention that, although by the virtue 
of the Central Limit Theorem, for any ^(77) with a given 
average v and dispersion D the distribution P(h,t) can 
be approximated by a Gaussian (||), the precision of this 
approximation is not sufficient to calculate averages of 
the type (W m {t)) = J e mh ^ P(h,t)dh. This integral is 
too sensitive to the precise shape of the distribution at 
the upper tail (or lower tail for m < 0). Indeed, the 
growth rate of ln(JT m (f)) equal to In dn ^(77) e^™, 
depends on the whole shape of 71(77) an d not only on its 
first and second moments v and D. 
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B. Multiplicative random walk in the continuous 
time approach 

The above multiplicative process is defined without 
ambiguity for discrete time. Straightforwardly taking 
the continuum limit causes problems. It might be use- 
ful to rewrite the equation of motion of a multiplica- 
tive random walk as a Stochastic Differential Equation 
(SDE) in continuous time. One should always keep in 
mind that a stochastic differential equation is nothing 
more than a convenient notation to describe a stochas- 
tic process in discrete time. At the n-th time step of 
discretized dynamics we define a new "continuous" time 
variable t as t — nAt. Here we introduced a rescal- 
ing factor At <C 1, which makes one step of underly- 
ing discrete dynamics an "infinitesimally small" incre- 
ment of the continuous time t. In the SDE approach one 
is limited to Gaussian distributed random variables, so 
we select a Gaussian distribution of 77(f) in our discrete 
dynamics. Since we want to approximate W(t) with a 
continuous function, the difference Wit + At) — W(t) 
after one step of discrete dynamics should be "infinites- 
imally" small. Therefore, we should select both the av- 
erage value and the dispersion of the Gaussian variable 
77(f) to scale as some power of At. It turns out to be 
the right choice to make them both scale linearly with 
At: 77(f) = vAt + Srj(t), where (Srj(t) 2 ) = DAt. Now one 
can write the discrete equation of motion for W(t) as 
W(t + At) = e vAt+Sr '(t) W ^ _ [i + v At + 67](t) + (vAt + 
5rj(t)) 2 /2]W(t) ~ W(t) + ((« + D/2)£t + 6r,(t))W(t), 
where we have dropped all terms smaller than linear in 
At <C 1. The SDE for W(t) can now be written as 

^p- = (v + §)W{t)+fj{t)W{t). (6) 

Here fj(t) = Srj(t)/At is a usual gaussian "continu- 
ous noise" with zero mean and correlations given by 
(fi(t)rj(t')) = D5(t — t'). We also assume the ab- 
sence of correlations between W(t) and fj(t). This as- 
sumption corresponds to selecting the Ito calculus over 
Stratonovich calculus. Both are just two formal ways of 
linking the polemic continuum limit and the well defined 
discrete version. 

The nontrivial part of this equation is an extra D/2 
term added to a deterministic growth rate of W(t). This 
term is not an artifact of our approach, but has a real 
physical meaning. Indeed, Eq. fl) can be solved for 
(W(t)) to give (W(t)) = (W(0))e$+ D M t , which is the 
right answer (see Eq. (|J)). Without this extra term 
we would be lead to the conclusion that for v = 0, i.e. 
(77) = 0, the average (not typical!) W(t) does not grow, 
which is wrong. 

The other way to get this extra term in the equation 
for W is to start with the well known Langevin equation 
of motion for h(t) — In W{t) describing a usual random 
walk with a drift: 

dh(t) _ 



where again (fj(t)) = 0, and (fj(t)ij(t')) = DS(t - f). 
To derive the equation of motion for W(t) = erW one 
has to do the change of variables as for usual partial 
differential equations. But in addition to this one has 
to add the "Ito term" || given by 4^ d Q%. 1 which is a 
formal prescription of Ito calculus. With this nontrivial 
correction one recovers the equation of motion (|^) . So in 
the language of SDE the difference between the typical 
(v) and the average (v + y) growth rates of W(t) in the 
multiplicative random walk is a direct consequence of the 
Ito term, appearing after the change of variables from h 
to W in the equation (^). 



C. Multiplicative random walk in the presence of a 
lower wall 

Much attention was devoted recently |2|,|| to the anal- 
ysis of the problem of "multiplicative random walk, re- 
pelled from zero" . In the economical context it was first 
introduced by Solomon et al. In a simplest case one 
has a multiplicative random walk with a Gaussian ran- 
dom variable 77, having a negative average v = (77) < 0, 
and the dispersion D. In other words the typical value of 
W{t) exponentially decreasing in time, while its average 
may or may not grow in time depending on the sign of 
v + D/2. In addition to this one has an "external force", 
pushing W(t) up and preventing it from falling below 
some predetermined constant. This external influence, 
which will be referred to as "lower wall" , should not sig- 
nificantly affect the dynamics for large W. One way to 
introduce a lower wall is to add an additional positive 
"source" term b into the RHS of Eq. (|^). The eqs. (^), 
and (Q) now become 



dW(t) 
dt 

dh(t) 
dt 



= {v + D/2)W(t)+ri(t)W(t)+b; (8) 
= v + r](t) + bcxp(-h). (9) 



As we see the lower wall in Eq. (||) has a property of 
being "short-ranged" in /z-space, i.e. its contribution to 
the SDE for h(t) can be neglected for large positive h. 
But for negative h the strength of the wall grows ex- 
ponentially and compensates the negative drift already 
at h = — hi(\v\/b). It is easy to convince oneself that 
this stochastic process eventually reaches a stationary 
state, characterized by a stationary probability distribu- 
tion P(h). In this stationary state the negative drift of 
h(t) is precisely balanced with diffusion combined with 
repulsion from the lower wall. 

In the literature on this subject one encounters many 
different realizations of the lower wall mechanism. For 
instance, one can introduce a more general term bW s 
into the RHS of M flq,|7|. In the equation for h this 



term becomes be^ s ^' 



rich for any S < 1 describes an 



dt 



(7) 



exponential lower wall qualitatively similar to (|^). In- 
deed, the "source" term in (||) is just a particular ex- 
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ample of this more general term with 5 = 0. On the 
other hand, the terms with b < and S > 1 describe 
an "upper wall" , preventing h from becoming to big. 
In this case, in order for a stationary state to exist one 
needs a positive drift of h pushing it up against the wall. 
In |Q] the lower wall is introduced "by hand" : in their 
simulations the authors simply do not allow h(t) to fall 
below a predetermined constant h m i n . In other words, 
h(t + 1) = min(/i(i) + r](t), h m i n ). Such "infinitely hard 
lower wall" can be described by a term bW s with very 
large negative S. Finally, Cont and Sornette || consider a 
case when the constant b itself can depend on time obey- 
ing a deterministic and/or stochastic dynamics. Except 
for pathological cases, where typical b(t) exponentially 
grows or decays in time, it does not qualitatively change 
the results, compared to a time-independent lower wall 

An interesting feature of a multiplicative random walk 
with a lower wall is that it generically gives rise to a power 
law tail in the distribution of W in the stationary state. 
We proceed by reviewing various derivations of this re- 
sult found in recent literature P,pj. As was explained 
above, the lower wall's only purpose is to make the pro- 
cess stationary by pushing the variable up whenever it 
becomes too small. The drift due to the wall can always 
be neglected for large enough h. In the region, where this 
approximation is justified one can write a Fokker-Planck 
equation, taking into account only the multiplicative part 
of the process, equivalent to diffusion with a drift in the 
/Vspace. The stationary solution of the Fokker-Planck 

equation should satisfy -d ^^ + -j d gffi^ = 0. It is 
easy to see that P(h) = A exp(2vh/ D) is indeed a so- 
lution. Since v < 0, it exponentially decays for positive 
h. The deviations from this form start to appear only 
at low h, where the presence of lower wall cannot be ne- 
glected. This "Boltzmann" tail of the distribution of h 
corresponds to a power law tail of distribution of W = e h : 
P(W) = AW~ 1+2v / D . The exponent of this power law 
tail 



t= l-2v/D= l + 2\v\/D 



(10) 



is greater than 1, so that there are no problems with 
normalization. In case of a lower wall of the form 
be~ h (sec Eq. (^)) one can write an analytic solu- 
tion of the Fokker-Planck equation valid for any h. 
It is the Boltzmann distribution with a Hamiltonian 
H(h) = be~ h — vh and temperature T = D/2, i.e. 
P(h) = Aexp[(-2be~ h /D + 2vh/D)], or P(W) = 
Aexp(-2b/DW)W~ 1+2v / D . The normalization con- 
stant A is given by A = (2b/ D)- 2v / D /T{-2v/ D). 

The Eq. ©, expressing the exponent of the power 
law tail of P(W) in terms of v and D, is valid only for 
the case of Gaussian distribution ^(77). Indeed, in its 
derivation we employed a stochastic differential equation 
approach, which is restricted to Gaussian noise. It is 
instructive to derive an equation, giving the value of r 
for a general ^(rj). It was first done by Kesten in @] 



and recently brought to the attention of physics com- 
munity in |3|]. Again, the formula holds for any mul- 
tiplicative process with a negative average drift ((7/) < 
0) and a lower wall, the effect of which can be ne- 
glected for large W. We assume that the process has 
already reached a stationary state, characterized by a 
stationary distribution P(W). For sufficiently large 
W, so that one can neglect the effect of the wall, the 
stationarity imposes the following functional equation 
on P{W): P(W) = J^n( V ) drj f+°° P(W) S(W - 
eW)dF' = dTjn^e-vPiWe-v). Assuming that 
the solution has a power law tail P(W) ~ W~ T one finds 

solution of 



1. In other words r is given by a 



= 1. 



(11) 



The obvious solution r = 1 should be rejected because 
the distribution function is not normalizable in this case. 
In short, we are looking for a solution with r > 1. Let us 
define A(r) = (e^ 7 "" 1 )). Since <2A(l)/dr = (77) < 0, but 
d 2 A(r)/dr 2 > one has at most one such a solution. In 
fact, if the distribution of p(rj) is not restricted to 77 < 0, 
for t — ► +00 one has A(r) — > +00 and the solution is 
guaranteed by the continuity of Kit). Only in the situ- 
ation when 77 is always negative, the region of large W 
is absolutely inaccessible, and no power law tail at large 
W is feasible. Using Eq. (fih, one can check that for a 
Gaussian distribution Eq. ( pd| ) predicts r = 1 — 2v/D in 
agreement with (10). 



D. Interpretation of W(i) as a fluctuating stock 
capital 

In what follows we will stick to the following "realiza- 
tion" of the random multiplicative process: we interpret 
W(t) as the capital (or wealth, hence the notation) that 
a single investor has in some stock. The price of the share 
of this stock p(t) undergoes a random multiplicative pro- 
cess p(t + 1) = e' ,< - t - ) p(i), and if the investor keeps a fixed 
number K of shares without selling or buying this stock, 
his capital W(t) — Kp(t) follows these price fluctuations. 
Later on we will consider models, where the investor at 
each time step will sell some stock and buy another. We 
assume that volumes of such transactions are sufficiently 
small, so that they have no influence on the market price 
fluctuations. Hence our assumption that rj(t) and W(t) 
are uncorrelated. 

The lesson one derives from the above properties of 
multiplicative random walk is that if the investor keeps 
all his money in just one stock it is the typical growth 
rate (77), he should be concerned about. In majority of 
realizations his capital grows at typical rate and he can- 
not directly take advantage of a bigger average growth 
rate ln(e' 7 ). There are situations when the typical growth 
rate is negative, i.e. the stock price is going down, while 



4 



the fluctuations are strong enough to make the average 
rate positive. The question we are going to address in 
this manuscript is how one can still exploit this average 
growth rate by investing and actively managing a port- 
folio composed of N stocks. 

III. ENSEMBLE OF TV STOCKS WITHOUT 
REDISTRIBUTION. 

The first problem we are going to consider is: what 
is the typical growth rate of the capital invested in an 
ensemble of N stocks if one is not allowed to sell one 
of them and reinvest the money into another. In the 
following we assume that the price Pi(t) of a share of 
each stock undergoes a multiplicative random walk, in- 
dependent of price fluctuations of other stocks. In other 
words, one time step logarithmic price increments rji{t) 
are uncorrelated not only at different times, but also for 
different stocks at a given time. The validity of this ap- 
proach for the real stock market lies beyond the scope of 
this work. For simplicity of final expression in this sec- 
tion we will restrict ourselves to the situation when rji for 
each of the stocks are Gaussian variables with zero mean 
((r/i) = v = 0) and the same dispersion D — (r)i). Ini- 
tially the capital is equally distributed between all stocks. 
We assume that the starting capital in each stock is equal 
to 1/N, so that the total capital is equal to 1. The typical 

value of the total capital (W to t{t))ty P = (J2iLi W / i(*))t yP 
will then grow in time. From the results of the pre- 
vious section one concludes that (Wj(£)) = e Dt / 2 and 
(W l (t) 2 ) - (W l {t)) 2 = e 2Dt - e Dt . One can safely re- 
place the sum of N variables with their average as long 
as (((VF i (<) 2 )-(VF i (t)) 2 )/iV) 1 / 2 < (Wi(t)). Therefore, at 
short times, when Dt <C IniV, one indeed enjoys the av- 
erage growth rate: (W{ t(i))typ = e Dt l 2 . At later times, 
however, the typical value of the capital starts to fall be- 
low the average value (i.e. average value over infinitely 
many realizations). To determine this slower growth of 
typical value quantitatively one has to approach the prob- 
lem from a different end. At late times the value of the 
total capital is mainly determined by the capital accu- 
mulated in the most successful stock, i.e. Wtot(i) — 
Wmax = maxj^jv Wi(t). The extremal statistics the- 
ory H readily gives the typical value of the W mM by 
requiring that 1/N = Prob(PF > W max ) = Prob(lnW / > 
In Wmax) ~ exp(— In 2 W max /2Dt) With exponential pre- 
cision one gets W max <~ e y ,2Dt ln N m Our approxima- 
tion that Wtot(i) — = maxi = i,jv W,(t) is good only if the 

(2) 

second maximal W (we denote it as Wmax(t)) is much 
smaller than the maximal one. Following the same argu- 
ments as before one concludes to find the typical value 
of w£L{t) one needs to sol ve Prob(VF > W^L) = 2/N. 
This results in W^L ~ e V 2£l *( ln N ~ 2 ). One easily con- 
firms that the approximation of the whole sum with its 
biggest element makes sense if Dt ^> IniV, which is 



a complementary condition to the "average" growth at 
small times. Therefore, we conclude that 

(WtotW)typ =e m ' 2 for t^lnN/D; (12) 
(W tot (t)) typ = e V2DtlnN for t^lnN/D, (13) 

Since growth proportional to \ft is slower than linear in 
t one concludes that no matter how big is your N your 
asymptotic growth of your total capital is still determined 
by the "typical" growth rate v = (77) (equal to zero in the 
case considered above) of a single stock. 

If one wants to exploit the "average" growth rate for 
a period of time T and then sell the stocks one needs to 
take an exponentially large ensemble of stocks N > e DT . 

IV. ENSEMBLE OF N STOCKS WITH 
REDISTRIBUTION. 

The case of "non-interacting" stocks, considered in the 
previous section, can be also called the case of a "lazy 
investor" . Indeed, initially the investor puts equal cap- 
ital in N stocks and leaves them as they are. He never 
sells or buys stocks. No wonder that very soon he can 
no longer expect to get an average rate of return on his 
investment and has to settle for smaller typical growth 
rate. Now we are going to consider the case of an ac- 
tive investor who after each time step redistributes his 
capital between stocks according to some simple rule. 
One may naively think that by selling unsuccessful stocks 
with small W% and reinvesting the money into successful 
stocks with large Wi one may do better. In reality the 
answer is precisely the opposite: one needs to sell some 
of the most successful stocks and reinvest the money into 
the least successful stocks. Selling only small number of 
shares of the most successful stocks (i.e. ones which are 
currently overpriced) and reinvesting this money into the 
least successful stocks (i.e. ones which are currently un- 
derpriced) makes a huge difference: ln Wi for underpriced 
stocks goes up significantly, while ln Wi for overpriced 
stocks does not go down as much. As we will show such 
a "charity" between stocks bootstraps the typical growth 
rate of the capital, so that ln(W to t(i))typ at all times has 
a growth rate bigger than a typical growth rate of a sin- 
gle stock. For large N this rate quickly approaches the 
average growth rate given by ln(e ,) ) (equal to D/2 for 
the Gaussian distribution of r\ with zero mean). This 
growth rate serves as a theoretical maximum of all pos- 
sible growth rates achievable by simple redistribution of 
funds. 



A. Problem with redistribution in the discrete time 
approach 

We start with the simplest strategy for redistribution 
of the capital. Under this strategy at each time step the 
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investor calcul ates t he current value of average capital 
per one stock W(t) — ^ E,=iAr^W- The capital is 
redistributed between the stocks according to the rule 
Wi — ► Wi — X(Wi — W). For positive A it means that 
"overpriced" stocks with Wi{t) > W(t) loose a fraction 
of their capita l in favor of the "underpriced" ones with 
Wi(t) < W(t). The extremal case of A = 1 corresponds 
to the equal redistribution of the capital after each time 
step. The stock price changes during the next discrete 
time interval. As a result the capital invested in each 
stock is multiplied by the random factor e 7 ^*'. The com- 
plete change of each stock's capital after one time step is 
given by: 



Wi(t + 1) = e"' (t) [(l - A)Wi(t) + XW(t)}. 



(14) 



One can recognize the above model can be inter- 
preted as the Directed Polymer model in N dimen- 
sions, with mean field (fully connected) interactions 
The role Laplacian is played by W{t) — Wi{t) 



l/N J2j=i n — It is convenient to intro- 

duce a new set of rescaled variables Si(t) — Wi{t)/W(t). 
The sum of Si is always equal to N , which sets a theo- 
retical cutoff equal to N to a value of individual s,-. One 
can rewrite Eqs. (|14|) in the following form: 



8i (t + 1) = ™® e^W[(l - \)a t (t) + A]; (15) 
W{t + l) = W(t) Z ^ l - 1 ' N ± — -. (16) 



As we will confirm later, the dynamics of W(t) can 
be approximated as a random multiplicative process, 
where the multiplication factor T(t) — Yli=i N e m ^((l — 
A)sj + X)/N has only small fluctuations around its av- 
erage value. We will indeed demonstrate that T(t) = 
(r) + 5r(t), where \ST(t)\ ~ N~ a / 2 . It means that 
for large TV to a good approximation one can disregard 
the fluctuations of W(t + 1)/W(t) while trying to solve 
Eq. (U). The average value of this ratio is easily 
calculated and is equal to (e n ) (one has to recall that 
n s i = N). In this approximation the equations 
of motion for Sj decouple and allow for exact solution. 
These mean-field equations are: 



e Vi(t) 



\{l-X) Sl {t) + X] 



(17) 



Similar equation of motions were recently studied by 
Cont et al. || and Solomon et al. Q and were shown 
to give rise to a stationary distribution of s having a 
power law tail for large s. One has to keep in mind that 
the definition of s in our problem introduces a natural 
cut off to this tail as s < N, so it is only for large N that 
one has a chance to see the effect of this power law or 
measure this power law numerically. 



The stationary distribution P(s) is conserved by dy- 
namics. Therefore, it should satisfy the following func- 
tional equation: 



P(s) = dr) n(rj)P( 



A 



1 - A 



)/R(tj), (18) 



where R(rj) = (1 — X)e v / (e v ) . Using this equation one 
can easily verify that indeed (s) = J sP(s)ds = 1, 
which is to be expected since X)i=i n Si = N. Assum- 
ing that P(s) has a power law tail of the form As~ T , 
and substituting it to the functional equation ( |l8| ) one 
gets the self consistency condition for the exponent r: 
J drj'K(ri)R(ri)' T ^ 1 = 1, or 



t]{t-1)\ 



1 



1-X 



(19) 



For a general distribution 77(77) this equation cannot be 
solved analytically. All one can deduce is that for a weak 
coupling Ac 1 the solution exists and is approximately 
given by t = 2. That means that for a weak coupling 
one always has P(s) ~ 1/s 2 ! For a case of Gaussian dis- 
tribution of 77 the analytic expression for r can be easily 
obtained from Eq. (0) and is given by 



21n(l-A) 
D 



(20) 



In Fig. 1 we present the results of simulations of the 
model with N = 10000. The measured power law expo- 
nent is in excellent agreement with the above theoretical 
prediction. 

Our ultimate goal is to determine W(i) t as a func- 
tion oft. The Eq. (Uq) states that at each time step W(t) 
is multiplied by T{t) = £i=i,iv e*W((l - X) Si + X)/N. 
One can show that T(t) at different time steps are un- 
correlated. One can a lso di sregard possible correlations 
between the value of W(t) and T(t) at the same time 
step. Then the behavior of W(t) is nothing else but a 
multiplicative random walk studied in Section 1. The 
typical value of W(t) grows as (W(t)) typ — e^ lnI ^ , while 
its average value grows as e tln ^ = e tln ^ e ^ = (e eta ) t . 

We will proceed by demonstrating that for any A > 
the typical and average growth rates of W(t) differ by 
0(N~ a ). For (lnT) one has the exact expression: 

(lnr) = ln(e") ^ln (\ + 1 X< [(1 - X)s t + A] J \ 

(21) 

where we introduced the notation Xi — e ,?i /(e , ') — 1. Ex- 
panding the second logarithm for large N, we get to lead- 
ing order: 



(lnr) ~ ln(e") - -L - X) Sl + A] 2 



N 



2N 2 



G 



where we used the fact that are uncorrelated at 

different i's. Therefore, (Xi(t)Xj(t')) = D5ijSt } t>, where 
D = (e 2, >) / (e^) 2 - 1. The fact that these variables are 
uncorrelated at different times proves that indeed W(t) 
undergoes a multiplicative random walk. The last step 
is to estimate Yli=i n s f • ^ ^° ^ n ' s we nee d to recall 
our results for the stationary distribution P(s). If the 
exponent r of the power law tail of this distribution is 
larger than 3, (s 2 ) is finite, J2i=i n s l = N(s 2 } and one 
immediately gets (lnT) = ln(r) — A/N, where A = [(1 — 
A) 2 (s 2 ) + A 2 + 2(1 - \)}D. In reality this is not hundred 
percent true. Indeed, expanding the logarithm in Eq. 

we stopped at the first order. In the presence of 
power law tails in P(s) the validity of this approximation 
is in doubt because the higher order terms involve the 
sum of powers s k with k > 2. For large enough k such 
powers are known to diverge as some power of N . It can 
be shown that for very large N they would dominate the 
scaling with respect to N. Such crossover was indeed 
observed in simulations. In Fig. 2 we present the results 
of the simulations of our model with A = 0.1, D = 0.1, 
which corresponds to r = 4.1. Indeed, we observe that for 
N — > oo, the difference between the average and typical 
growth rates of the total capital, v avg — v t yp(N) = ln(r) — 
(lnT), approaches zero. This approach starts as A/N a 
with a = 1, but at larger N a deviation towards smaller 
a can be noticed. 

For 2 < t < 3 the second moment of s diverges. 
This means that one should be more careful in estimat- 
ing J2i=i n s l- The apparent divergence of the integral 
J s 2 P(s) ds should not be taken too seriously, since we 
are dealing with a finite sample of variables s restricted 
by ^2 Si — N. Even in the worst case if only one 
is nonzero (and equal to N by normalization) the sum 



= 1,N ■ 



= N 2 



In all situations when the integral 



J s k P(s) ds diverges the sum of a finite sample is dom- 
inated by the largest element. One can estimate this 



largest s by requiring Prob(s > s n 



l/N. 



Therefore, the typical value of the largest Sj is given by 
Smax ~ iV 1 /^ 1 '. Since r > 2 this value is always less 
then N - the maximal possible s. Then X)i=i n s f — 



A, v, and D of a discrete model very small, keeping their 
ratio fixed. 

In the limit At <§; 1 the Eq. ( |l4| ) becomes a stochastic 
differential equation 

d t Wi(t) = \(W-Wi) + (v + D/2) Wi + Wifji(t). (22) 

Here as in Section 1 we introduced the continuous-time 
stochastic force fji(t) — r)i(t)/At — v, and used e Vi = 
l + m + r 1 2 /2 + ...~l + f H At+(v + D/2)At + 0(At 3 / 2 ). 
It is important to point out here that such a continuous 
time formulation is only meaningful if rji(t) is a Gaus- 
sian noise. Only in this case Eq. (|2^ ) can be regarded 
as a Langevin equation Q . Usually the assumption of a 
Gaussian noise is motivated by the fact that for a con- 
tinuous time process, the stochastic force fjidt acting on 
a small interval At can be thought of a sum of infinitely 
many infinitesimal contributions. The central limit the- 
orem then ensures that rjiAt is Gaussian. For processes 
with additive noise, this assumption is reasonable also 
for discrete time processes. For multiplicative processes 
the deviations from the central limit theorem becomes of 
concern since the tails of the distributions are probed by 
the process. Therefore, we shall assume in this section, 
that fji is Gaussian. 

Under this assumption, we shall be able to derive the 
full probability distribution of the W% in the limit N — > 
oo. It is again convenient to use the variables s$(t) = 
Wi(t)/W(t). Using Ito calculus, one readily finds 



d t s l = A(l - 



D 



Si 



+ Sifji - sfj Si (23) 



where we used the notation f = jf Ylj fj ■ Note that 
J2i Si = N and, consistently, J2i &t s i — 0- We shall 
adopt a self consistent mean field approach, valid in the 
N — > oo limit, in which we substitute averages over i with 
statistical averages: /= (/). Within this approximation, 
the term sfj = (sfj) = can be neglected. If we introduce 



2A 2 . 2x _ 2A 2^ 

(s ) = s 2 

D N D N 



(24) 



N 2 /( 7 



as a constant to be determined later self-consistently, 
Eq. (B3|) becomes an equation for Sj only, which does not 



comes (lnT) = ln(T)-A'/N 



Now the expression for the (lnT) be- invol ^ . for • ^ l explicitly . We know S that, for a 



with A' ~ (1— A) D. Langevin equation of the form 



B. Problem with redistribution in continuous time 
approach 

Similar results can be obtained in the continuous time 
limit of Eq. (|l4]). In order to derive the stochastic par- 
tial differential equation corresponding to Eq. (^) we 
assume that time is discretized t — nAt in units At and 
we take A = A c Ai, v — v c At, and D = D c At. In all our 
future formulas we drop the subscript c in A c , v c , and D c 
of continuous model. However, one should keep in mind 
that we recover continuous limit by making parameters 



d t s = 



> dV(s) 
ds 



ST], 



the associated Fokker Planck equation yields the asym 
totic distribution P(s) ~ e 
into this form, we find V(s) 
which 



2V(s)/d_ Recasting Eq. Q23|) 
= X/s + ^^- + Ds/N, from 



P(s) = Afexp 



2A 
Ds 



N 



(25) 



Note the emergence of a power law behavior in P{s), 
which is however cut off by the second term in the expo- 
nential. This is physically meaningful, since s < N must 



7 



hold, with s — N occurring when the whole capital NW 
is invested in a single stock. The value r of the power 
law decay is determined self-consistently from Eq. (H ) 
performing the average on the distribution in Eq. ( |25| ). 
A further requirement which our approach imposes on 
P{s) is that s = (s) = 1. It is not possible to compute 
exactly these averages, however, it is possible to perform 
a large N expansion. Indeed if we set 



Z(p) 



ds exp 



2A 
~Ds 



then 

32 



and 



clearly, (s) = -d fl \nZ(fj,)\ fl=2 /N 
9^1nZ(^)| A1=2 /Af + (s) 2 . Therefore, evaluating the small 
fi expansion of Z(n) we can compute the first two mo- 
ments of s and impose self-consistency. However Z(p) 
has a non-analytic expansion around fi = 0, since deriva- 
tives d™Z(fj,) diverge at \i — for n > r. For A > D/2, 
the first two derivatives exist. The equation (s) = 1 then 
allows us to compute r together with its leading correc- 
tion: 



„ 2A 
r*2 + - 



X{{D/2) 2 - XD/2 + X 2 ) 4 
(D/2) 2 (D/2 - A) N : 



forD/2 < A 
(26) 



The equation ( p4| ) then turns out to be automatically 
satisfied, which is a reassuring check of self-consistency. 
Note that Eq. ( |20| ) derived previously, exactly reduces 
to Eq. ( |26| ) with A <C 1. For r < 3 the second derivative 
of Z(/i) does not exist at /i = 0. The second term in Eq. 
(p6|) changes, but the leading term remains the same: 



2A 2X/D + 1 
2 + — --_ 



D r(2A/D + l) 



dx- 



1 



r 2+2A/D 



4A 



The average growth rate of the capital NW is obtained 
summing Eq. (E2T) over i and dividing by N: 



d t W{t) = (v + D/2 + sf)) W. 
The solution to this equation 



W{t) = M/(0)exp 



(v + D/2)t+ / sfj(t')dt' 
Jo 

^F(0) e (" +fl / 2 ' ( 



implies that the growth rate of the average is, to leading 
order in N, equal to the growth rate of the average v + 
D/2. 



C. Parallels to directed polymers in random media 

In conclusion we would like to point out that 
the stochastic differential equation (|22]) has a finite- 
dimensional analogue, which was much studied over the 
past decade. Indeed, the term X(W — Wj) is nothing 



else but a fully connected (infinite dimensional) variant 
of discrete Laplacian. In finite dimensions this term be- 
comes XAWi = X(J2 nn W R n/2d - Wi). In the spatial 
continuous limit the Eq. (|22J) becomes 



d t W(x, t) = XAW{x, t) + (v + D/2) W(x, t) 

+ T](x,t)W( X ,t), 



(27) 



which can be easily recognized as the equation for the 
par tition function of directed polymer in random media 
|10f . The change of variables h = In W maps this equa- 
tion to the so-called KPZ equation [fTT|| : 

d t h(x, t) = X(Ah(x, t) + \Vh(x, t)\ 2 ) + vh(x, t) + rj(x, t) 

(28) 

In our infinite-dimensional (fully connected) model we 
found that P(W) has a power law behavior for large W . 
In finite dimensions, at least below the upper critical di- 
mension d c (whose very existence is still under debate), 
this seems not to be the case. Indeed numerical simula- 
tions show that, at least up to d — 3 + 1 jl2| the distribu- 
tion of h = hiW has not a pure exponential, but rather 
stretched exponential behavior. We conjecture that the 
power law behavior of P(W) in the model studied in this 
manuscript is an artifact of the peculiar long range inter- 
action, where each site is coupled to any other site. 

The work at Brookhaven National Laboratory was sup- 
ported by the U.S. Department of Energy Division of 
Material Science, under contract DE-AC02-76CH00016. 
S.M. and thanks the Institut de Physique Theorique for 
the hospitality during the visit, when this work was ini- 
tiated. This work is supported in part by the Swiss Na- 
tional Foundation through the Grant 20-46918.96. 



[1] 

[2] 

[3] 
[4] 

[5] 

[6] 

[7] 

[8] 

[9] 
[10] 



R.C. Merton, Continuous-Time Finance, Blackwell, 
Cambridge (1990). 

Levy and S. Solomon, Int. J. Mod. Phys. C7, 595 (1996), 
ibid, 745 (1996). 
R. Cont and D. Sornette, 

B. Gnedenko "Theory 
Moscow) 

C. W. Gardiner, Handbook of Stochastic Methods, second 
edition, p. 124 ( Springer- Verlag 1985). 

Y. Tu, G. Grinstein and M.A. Mufioz, Phys. Rev. Lett. 
76, 4376 (1996). 



J. Physique I, 7, 431 (1997). 
of Probability" (Mir 1978, 



M.A. Murioz and T. Hwa, lcond-mat/9702217] preprint. 
H. Kesten, Random difference equations and renewal the- 
ory for products of random matrices, Acta Math. 131, 
207 (1973). 

J. Galambos, The Asymptotic Theory of Extreme Order 

Statistics (John Wiley & Sons, New York 1978). 

T. Halpin-Healy and Y-C. Zhang, Phys. Rep. 254, 

(1995). 



8 



[11] M. Kardar, G. Parisi, Y.-C. Zhang, Phys. Rev. Lett. 56, 
889 (1986). 

[12] L. M. Kim, A. J. Bray and M. A. Moore, Phys. Rev. A 
44, 2345 (1991). 

FIG. 1. The distribution of capital fractions Si = Wi/W 
for A = 0.25, 0.5, and gaussian tt(t]) with D — 2, and v = in 
a system of size N = 10000. The solid lines are the theoretical 
predictions (20) for a power law exponent r of the tail of this 
distribution. 

FIG. 2. The difference between the average growth rate of 
the capital v avg — D/2 and its typical growth rate v typ (N) as 
a function of the number of assets N. The parameters of the 
model are A = 0.1, D = 0.1, v = 0. The solid line indicates the 
theoretical prediction A/N. The crossover towards smaller a 
is clearly seen for large N. 
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Fig. 2 




system size N 



